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We investigate anharmonically coupled high- and low-frequency excita- 
"^ . tions in a planar lattice of adsorbed molecules interacting with phonons of a 

crystal. Dephasing of high-frequency local vibrations by low-frequency res- 
onance modes is described in terms of temperature Green's function. The 
equations obtained are solved, first, with a small ratio of the anharmonic cou- 
^ . pling coefficient for high- and low-frequency modes to the resonance width, 

o. 

^ ■ and second, in the low-temperature limit. High-frequency spectral line posi- 

tions and widths depend on dispersion laws and resonance mode lifetimes. It is 
^L . shown that lateral interactions of low-frequency modes of adsorbed molecules 

\^ . can lead to a significant narrowing of high-frequency spectral lines, which is 

o, 

On , consistent with experimental data. 

I ' A considerable progress in the understanding of the broadening mechanism for a high- 

Tj ■ 

Q I frequency local vibration spectral line of an impurity molecule has been achieved due to 

O ■ 

studies on the self-scattering (dephasing) of this vibration arising from a biquadratic anhar- 

•T-H , 

/\ ' monic coupling between high- and low-frequency molecular modes [1-3]. At a sufficiently 

high frequency of the local vibration, other anharmonic coupling types do not contribute 
to the spectral line broadening [3-5] with the exception of additionally renormalizing the 
biquadratic anharmonicity coefficient [6]. In this context, of significance is that the low- 
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frequency mode acquires a resonance nature due to its interaction with the quasicontinuous 
spectrum of crystal lattice phonons. In the treatment of the self-scattering effect, the de- 
phasing model [4,7,8] proved fruitful, as it afforded a number of exact solutions [9-12] and 
hence a description of spectral line shapes for local vibrations. 

If the amount of adsorbed molecules suffices to form a monolayer, the dephasing model 
should take into account collectivization of vibrational molecular modes [13-16]. As a rule, 
high- and low-frequency lattice excitations are respectively related to stretch and deforma- 
tion (translational or orientational) molecular vibrations. 

In this paper, a dephasing model for local vibrations is investigated for a system of 
anharmonically coupled high- and low-frequency molecular vibrations, with crystal phonons 
included. If translational symmetry is characteristic of the system considered, the low- 
frequency contribution to its Hamiltonian can be introduced in a diagonal form with respect 
to the two-dimensional wave vector K belonging to the first Brillouin zone of a planar 
lattice [17]. Then the total Hamiltonian is written as a sum of harmonic and anharmonic 
contributions: 

Htot = H + Ha-, H = ^ i^K, (1) 

K 

i^K = ^f^Ka^OK + hu-Kb^bK + hJ2^K,,yPK,uPK,'y + ^H(^K,i.&K/5K,z. + H.C.), (2) 

u V 

fl^ ^ 

Ha = "77 zJ flK'^K'&K"^K-K'+K", (3) 

^^ KK'K" 

where the terms in the expression (2) respectively represent high- and low-frequency vi- 
brations of a molecular lattice, crystal phonons, and the harmonic interactions between 
phonons and low-frequency modes. In the above relation, quantum states of phonons of 
the admolecular lattice are characterized by the surface-parallel wave vector K, whereas the 
quantum numbers of substrate phonons are indicated by the couple of indices K and u. The 
latter accounts for the polarization of a quasi-particle and its motion in the surface-normal 
direction; it implicitly reflects both the atomic arrangement in the crystal unit cell and the 



position of the admolecular lattice, as a whole, relative to the crystal surface. The utility of 
introduction of the complex index v was substantiated previously [17]. 

As a rule, the density of states for molecular lattice vibrations is negligible as compared 
to that for crystal phonons. Therefore, the K-mode of a molecular lattice interacts with 
the crystal phonons specified by the same wave vector K. Besides, the low-frequency col- 
lective mode cuk of adsorbed molecules can be considered as a resonance vibration with the 
renormalized frequency cI^k and the inverse lifetime r^K ('^k) [12,17]. 

The response of the system concerned to an external electromagnetic field is conveniently 
described in terms of double-time Green's function (GF) which can be introduced in a 
variety of representations [18-22]. In what follows we will involve the representation in 
Matsubara's frequency space [18] which is accepted in the theory of anharmonic crystals [23] 
and provides a number of exact solutions in the case of a single adsorbed molecule [10,11]. 
In this approach, the spectral line shape for high-frequency vibrations can be determined 
as follows: 
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and A is an arbitrary operator function. 

Expanding the right side of Eq. (6) in powers of the anharmonic couphng coefficient 
7, the temperature GF (6) is expressed in terms of the GF's of high- and low-frequency 
vibrations of a planar lattice in the harmonic approximation: 

= {[n{Qj^) + l]e{r) + niQj^)e{-r)}exp{-hQj^r), (9) 

9^K\r) = (f6K(r)6+(0))^ 
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where 

n(cu)= [e'^^^-lj"' (11) 



is the is Bose-factor and 



^K{io) = -, z —2 : — (12) 



is the resonance spectral function of low-frequency vibrations [12,24] 
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The real and imaginary parts of the pole in expression (12) define the renormalized 
frequency ujk and the inverse lifetime r^K = ^Trfjici^K.) of the resonance vibration. 

Let us take advantage of the inequality (3hQ-K » 1 permitting neglect of the terms of 
the order exp{—f3hQ-K)- Then the trace taken over high- and low-frequency modes, Sp{...}, 
is reduced to that for low-frequency modes, with all high-frequency vibrations considered 
only for the ground state. The resulting temperature GF (6) takes the form: 

G{ao4, r) = (ao(r)^(r)5+(0))^, P > r > 0. (14) 



Thus, provided the inequahty f3hQti >> 1 is vahd, the GF defined above can be written 
in the form involving no high-frequency mode operators: 

G{aoa^, t) = exp (/il^o^) (soo(r))^, (15) 

n=0 { •'q Ki...K„_i 

^K(r) = ^exp [h{ny, - l^K')] E^K"-K(r)&K"-K'(^)- (17) 

Representing average value (15) as 

'Mr))^ = e^(^) (18) 

and differentiating the left and the right sides of expression (18) with respect to r, we arrive 
at: 



!^-^ft,^i^^^<ll!!^. (19) 
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Eq. (19) allows the line shape to be determined in the approximation of the small an- 
harmonic coefficient 7, as well as in the low-temperature limit. 
Importantly, the three-time GF 



, „^ r6K(r)6i,,(r')sK-K',o(r"),„ 
gKK'ir, t', t") = ^ J — ^ ^ (20) 




satisfies the Dyson-like equation with the assumption of non-interacting high-frequency 

collective vibrations {Qk = ^o)'- 

9KK'{r, t', r") = 4°V - t')5kk' + ^ E / ^^i^kV " n)^K,K'(ri, r', r"). (21) 

Ki 

(the resonance GF giQ (t) in the harmonic approximation is determined by the expres- 
sion (10)). Equation (21) takes into account intermolecular interactions of low-frequency 
resonance modes and thus generalizes the corresponding equations derived in papers [10,11]. 



To invoke the perturbation theory for small anharmonic coupling coefficient, we use the 
Wick theorem for the coupling of the creation and annihilation operators of low-frequency 
modes in expression (19). Retaining the terms of the orders 7 and 7^ we are led to the 
following expressions for the shift AQ and the width 2r of the high-frequency vibration 

spectral line: 

00 
^^ = ^ E / dujn{uj)^K{oo), (22) 
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where n{uj) and Kk(^) are respectively determined by relations (11) and (12). It should 
be noted that expressions (22) and (23) describe a single Lorentz-like spectral line of local 
vibrations, with its position and width dictated by dispersion laws and lifetimes for resonance 
low-frequency modes. 

The analytical description of high-frequency line shapes becomes possible in the low- 
temperature limit, i.e., at n{ujj^) ^ exp{— jShuic) << 1, which represents an experimentally 
important case. In this situation, the Wick coupling for the operators of low-frequency 
modes in expression (19) involves only the terms in which the annihilation operator is to the 
left of the creation operator in all but one operator pair. Then Eq. (19) can be written as: 
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On integrating Eq. (24) over the variables Ti, r2, . . . , t„, we obtain: 
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Then expression (26) in the Markov approximation [hriji_{uj)T, h'-fr » 1] can be reduced as 
follows: 

dr N ^l- {^/N) Ek' {"kk' + iTr [r7K(^K + %'(^k']}"' ' 

where 

ttKK' = ^0 — ^K-K' + ^K — ^K'- (28) 

As seen, the spectral line of high-frequency local vibrations is of the Lorentz-like shape: 

^(^) = --Im ^--, W = -^, (29) 

n uj — ilo — W dr 

with the shift Afi = ReW and the width 21 = -2lmW. 

To estimate the lateral interaction effect on basic spectral characteristics of high- 
frequency local vibrations, assume that the density of states for the low-frequency band 
is described by the step function: 

r Acj"\ uJo<uj <uJo + Au;, 
P^ (^) = (30) 

'^ 0, uj < LdQ or uj > uJo + Au, 

with the step width Au, which is normally much greater than the width of the high- 
frequency vibration band. In addition, suppose resonance widths ^k(^) to be much the 
same on the interval Ac<j [^k(^) ~ VK-i^o) = v]- Then, relation (27) takes the form: 



W = ^e-^-o/'^^^ / dxe-^^'^^"" 1 + T^ In "" ^ '"^ ~ . (31) 

AVL J \ An x + iT] J ^ ^ 

Assuming the sufficiently small ratio Au/rj and allowing for the fact that the quantities 
AQ and 2r essentially differ from zero only when /c^T approaches Uq (i.e., provided the 
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inequalities Ac<j << r] << ksT/h < ujq hold), we obtain, accurate to terms of the orders 

Au/T and {Au/r])^: 
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At Ac(j = 0, relations (32) and (33) are reduced to the reported results [8] for an ad- 
sorbed molecule. Lateral interactions of adsorbed molecules give rise to the low-frequency 
molecular band with the non-zero parameter Auo. The integration over this band leads to 
the decreased shift and width of the high-frequency spectral line. Since the parameters 
hAuj/kBT and [Auj/i])'^ in the temperature range of interest can be of the same order of 
magnitude, the contribution {Auo/rjY to the value AVt can be only positive, whereas the 
sign of this contribution to the value 2r depends on the 77 to 7 ratio. At 77/7 > 1.094 
(which is consistent with the realistic values of 77 and 7), the contribution {Auo/rjY causes 
the additional narrowing of the spectral line. 

Table 1 shows temperature dependences of the high-frequency line shift and width, with 
the parameter values corresponding to isotopically diluted ^^C^^02 molecules in the ^2c^^02 
monolayer on the NaCl(lOO) surface [25]. The data outlined are indicative of changes in 
spectral line shifts and widths when passing from the band width Auo = 5 cm^^ to Auj = 0. 
These dependences evidently demonstrate the narrowing of the spectral line, in accordance 
with experimental results [25]. In the framework of the conventional exchange dephasing 
model, i.e., with no regard to lateral interactions of low-frequency modes, an attempted 
description of the temperature dependences for the above-mentioned '^^C^^02 molecular 
ensemble resulted in overestimated line widths, whereas theoretically calculated line shifts 
agreed well with experimental values. Stretch vibration frequencies for ^^C^^02 and ^2c^^02 
molecules differ by 60 cm~^ and their coupling can thus be neglected. In contrast, low- 
frequency vibrations prove to be essentially coupled: As far as translational vibrations are 
concerned, the mass difference between carbon isotopes is slight compared to the total 

8 



molecular mass, whereas for orientational vibrations, the mass of the central carbon atom has 
only a slight effect on the molecular moment of inertia. With the coupling of low-frequency 
vibrations included, the model suggested provides an accordance between calculated and 
observed spectral line widths. 
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Table 1. Temperature dependences of the spectral line shift (a) and width (b) for 
molecular local vibrations, with intermolecular interactions neglected {Auj = 0) and in- 
cluded {Auj = 5 cm~^) for the following parameter values: 7 = 0.7 cm~^, rj = 1.2 cm~^, 
ujo = 41 cm"-"^. 

(a) 



T(K) 


AQ (cm-i) 
at Acu = 


AQ (cm-i) 
at Auj = 5 cm^^ 


20 


0.027223 


0.029096 


30 


0.072878 


0.082131 


40 


0.119241 


0.138051 


50 


0.160221 


0.188547 


60 


0.195096 


0.232113 


70 


0.224566 


0.269276 


80 


0.249553 


0.301008 



(b) 



T(K) 


2r (cm-i) 
at Acu = 


2r (cm-i) 
at Auj = 5 cm^^ 


20 


0.031761 


0.016631 


30 


0.085025 


0.046801 


40 


0.139113 


0.078537 


50 


0.186925 


0.107158 


60 


0.227612 


0.131831 


70 


0.261993 


0.152864 


80 


0.291145 


0.170817 
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